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r^ ' Abstract. 

$Ph' In this sequel paper we explore how macroscopic quantum phenomena can be measured or 

understood from the behavior of quantum correlations which exist in a quantum system of many 
^ ' particles or components and how the interaction strengths change with energy or scale, under 

2 I ordinary situations and when the system is near its critical point. We use the nPI (master) 

effective action related to the Boltzmann-BBGKY / Schwinger-Dyson hierarchy of equations as 
a tool for systemizing the contributions of higher order correlation functions to the dynamics 
of lower order correlation functions. Together with the large A'^ expansion discussed in our 
first paper [I] we explore 1) the conditions whereby an H-theorem is obtained, which can be 
^ ' viewed as a signifier of the emergence of macroscopic behavior in the system. We give two 

OO , more examples from past work: 2) the nonequilibrium dynamics of A^ atoms in an optical 

y^ • lattice under the large Af (field components), 2PI and second order perturbative expansions, 

illustrating how N and Af enter in these three aspects of quantum correlations, coherence and 
coupling strength. 3) the behavior of an interacting quantum system near its critical point, 
C^^ ' the effects of quantum and thermal fluctuations and the conditions under which the system 

^3 , manifests infrared dimensional reduction. We also discuss how the effective field theory concept 

bears on macroscopic quantum phenomena: the running of the coupling parameters with energy 
or scale imparts a dynamical-dependent and an interaction-sensitive definition of 'macroscopia'. 



1. Key points: Correlation, Interaction Strength, Effective Scales and Dimensions 

We continue our investigations begun in an earlier paper [Ij (MQPl) into the key features 
of macroscopic quantum phenomena (MQP), seeking to identify and assemble the necessary 
ingredients toward the construction of a viable theoretical framework. Not only does this new 
field bear on the foundational issues of both quantum and statistical mechanics it also sees a 
widening range of manifestations and applications. (For background please refer to references in 
[1]). We want to examine the conditions or criteria whereby a macroscopic quantum system may 
take on classical attributes, and, more interestingly, that it keeps some of its quantum features. 
Our first paper focused on the large N approximation, highlighting the subtle yet important 
difference between results obtained from a leading-order large-A^ approximation, namely, the 
mean field theory which is quantum in nature, from a corresponding classical theory. Here we 
focus on the key features of correlations and couplings in a quantum system of many particles or 
components, especially its behavior near the critical point, to investigate which factors contribute 



to, or determine a system as being macroscopic, and under what conditions would quantum 
features persist or dominate in that hmit. 

Before delving into these main subjects, it may be useful for the purpose of orientation to 
refresh what is already known in quantum statistical mechanics of a) how the thermodynamic 
properties of a quantum system depend on N, the number of particles present, and b) how 
qualitative changes arise in a many-body quantum system, such as the formation of a Bose- 
Einstein condensate (BEC), one of the most commonly cited examples (the other being quantized 
fluxes in Josephson junctions superconductivity) of MQP. We want to point out the qualitative 
differences in the meanings of 'quantumness' in these familiar cases before we launch our 
investigations into the other quantum features we are more interested in, namely, quantum 
correlations, fluctuations, coherence and entanglement, and how they may show up or even 
persist at a macroscopic scale. 

1.1. N dependence in quantum statistical mechanics 

Note from the start that one of the two basic premises of quantum statistical mechanics (QSM) 
[2], namely, the assumption of random phase approximation, precludes quantum coherence 
considerations. All descriptions are in terms of probabilities, not amplitudes. The other basic 
premise, that the system has equal a priori probability to be found in any of its accessible 
states, refers to the equilibrium condition. Note, however, that decoherence and quantum to 
classical transition [3l H] are fundamentally nonequilibrium processes. 

The 'quantum' in QSM refers to effects due to spin-statistics, in the difference between 
bosons and fermions, and distinguishability: different combinatorics in distinguishable 
(classical feature) versus identical (quantum feature) particles, which underlie the Gibbs paradox 
[5]. Thus the answer to a) above is simple: the dependence of the partition function Z on 
A'^ is explicit, and the thermodynamic limit is well known. With the correct counting of 
distinguishable versus indistinguishable particles included in the partition function, the difference 
between quantum and classical features in this restricted sense is also well accounted for. 

QSM describes finite temperature equilibrium processes for a quantum (in the above 
restricted sense) system in or near equilibrium, such as in the ordinary (thermal) critical 
phenomena mediated by energy fluctuations in a canonical ensemble (temperature), or numbers 
/ species fluctuations in a grand canonical ensemble (chemical potential) description, but cannot 
address issues related to quantum coherence and entanglement, the central issues in quantum 
information. To the extent entanglement is a useful marker for quantum phase transition QSM 
cannot address issues of quantum noise or fluctuation phenomena. 

BEC is the well-known process where one can pinpoint the number Nc of a system of N 
bosonic atoms undergoing a transition where the salient features of the system undergo marked 
qualitative changes on the large scale, exemplifying MQP. There, the critical number Nc and 
temperature Tc which define the phase transition also depend on the dimensionality. (E.g., a 
two-dimensional Bose gas undergoes transition into vortex-antivortex pairs - see e.g., [6]). BEC 
is certainly a quantum entity, carrying macroscopic quantum coherence described by the A^-atom 
wave function as many exquisite experiments carried out in the last decade have witnessed. Yet 
its pivotal feature, its 'quantumness', originates from the spin-statistics properties (bosonic vs 
fermionic condensate) and the critical phenomenon which describes the formation of a condensate 
is a classical, not a quantum phase transition. This addresses query b) above. 

Making the distinction between 'quantumness' arising from particle spin-statistics as in 
quantum statistical mechanics and that of large scale quantum coherence as in BEC offers a way 
to discern what would be called MQP in the modern sense of the word: usually it is the latter, 
not the former, which is the deciding factor. To name two examples: 1) Can the behavior of a 
collection of quantum gas be regarded as MQP? In solid state physics one treats the conduction 
electrons in a metal as a degenerate Fermi gas because the room temperature is much lower 



than the Fermi temperature of the metal. The electron's de Broglie thermal wavelength is much 
larger than the atomic spacing. But one would not regard the thermal properties of metals 
as MQP. Similar situations occur in the white dwarfs and neutron stars where the degenerate 
pressure of electrons or neutrons respectively balances the gravitation attraction to keep the 
stars in relative stable configurations. One does not refer to the existence of white dwarfs or 
neutron stars as examples of MQP. 2) Can one view collective excitations as MQP? Phonons 
are the quantized long wavelength modes of lattice vibrations. Like waves all such collective 
variables involve a large number of atoms and can span a large spatial extent. But one would 
not call these collective phenomena as MQP. However for a large population of bosonic gas 
accumulating in the ground state at an extremely low temperature forming a Bose-Einstein 
condensate(BEC), we do regard it as a MQP. The crucial difference is the large scale quantum 
coherence established in such systems. Similarly, for superconductors: the Cooper pairs existing 
at large spatial separation (compared to the extent of electron wave functions) are quasiparticles 
with full quantum coherence Q| 

1.2. Scaling and Interaction, RG running and Criticality 

In parallel to explicating the (rather restrictive) meaning of 'quantumness' used widely in QSM 
mentioned above here we wish to explore the meaning of 'macroscopia' in like spirit. Macroscopia 
in the most direct way conjures large numbers: many particles or components and/or large sizes 
or scales. Paper I focused on the 'many' aspect, here we focus on the 'scale' aspect. If one goes 
by the scale of a quantum system alone without other considerations such as the interactions (or 
coupling strength) amongst its underlying constituents then the quantum behavior of a macro 
object such as the Universe would be similar to that of a micro one, such as the internal degrees 
of freedom of an atom. In fact, both obey a harmonic oscillator equation of motion, albeit the 
former (for the wave function of the universe) has a negative spring constant. The scale is usually 
measured by the inverse mass of the relevant process. The physically relevant scale (or mass) 
is different from the bare scale (or mass) by a ratio given by a renormalization constant which 
captures the relevant interactions involved. Interaction is measured by the coupling constants 
which vary ('run') with energy according to the renormalization group (RG) equations. The 
infrared behavior of an interacting quantum system at or near its critical point such as the 
universality class it belongs to, is of special significance. We will review briefly these familiar 
concepts below as the notion of 'macroscopic' for quantum systems near criticality is determined 
by its infrared behavior and the interaction strength of its constituents. 

1.2.1. Scales, Interaction Strength and Effective Theories Let us examine the calculation of 
the energy levels of a hydrogen atom as an example. In a standard textbook treatment one uses 
the Schroedinger equation for an electron moving in the static Coulomb field of the proton. To 
a good approximation, the only properties of a proton which are relevant to this problem are 
its mass and charge. The knowledge of the quark structure of the proton is not necessary to 
compute the energy levels of the hydrogen atom. 

The required knowledge of the proton depends on how accurate one asks for the energy levels. 
A more detailed calculation, taking hyperfine splitting into account, for instance, requires the 
knowledge of the spin of the proton and the value of the magnetic dipole moment. An even 
more accurate calculation requires the knowledge of the proton charge radius and details of the 
proton structure. 

^ Another more practical angle or functional criterion towards this demarkation is to ask whether one could 
use such a system to perform quantum information processing, where quantum coherence and entanglement are 
essential. Cold atoms in an optical lattice, atom assembly, as well as superconducting flux qubits are the well- 
known viable examples. Note, however, there are proposed QIP schemes such as the NMR based ones, which do 
not invoke quantum entanglement directly, yet could achieve speed-ups over classical computation. 



The typical length scale characteristic of the hydrogen atom is the Bohr radius tq = l/{m(,a) 
(in units where fi = c = \ ). The typical momentum scale is 1/ro = meO, the typical energy 
scale characteristic of the hydrogen atom is ~ meOp' , and the typical time scale is l/(mea^). We 
can get a quantitative estimate of the errors caused by the neglected interactions. 

In an effective theory description the relevant interactions also depend on the question being 
asked. For the hydrogen atom the energy spectrum can be computed to the accuracy (mea/M^)'^ 
while ignoring the weak interactions. But if one is interested in atomic parity violation, the weak 
interactions are the leading contributions because the electromagnetic and strong interactions 
conserve parity. The effect of atomic parity violation will be very small because the weak scale 
is much larger than the atomic scale. 

As one approaches shorter and shorter distances (or higher and higher energies), one needs 
to invoke the relevant physical processes at that scale. The high energy processes affecting 
low energy physics appear in physical parameters measured at the low energy regime. This is 
captured nicely by effective field theories, the essence of which we will discuss in a later section. 

1.2.2. RG running and Criticality In an interacting quantum system the variation of 
interaction strength with energy is described by the RG equations which are derived from the 
renormalization of the coupling constants (for ultraviolet divergence) . Assuming that the system 
scales homogeneously this enables one to obtain the infrared properties of a quantum system 
from its ultraviolet behavior. At the critical point the correlation functions of such a system 
usually diverge: the system 'feels', so to speak, all the way to infinity. Let us consider how this 
interesting physical situation bears on macroscopic quantum phenomena. 

Critical phenomena involves many different length scales. Consider a physical system going 
through a continuous phase transition under some temperature change: the correlation length 
diverges as the temperature approaches the critical temperature, and fluctuations at all length 
scales are involved. When the temperature lowers across the critical point, the system's order 
parameter will develop a nonzero value, where the relevant degrees of freedom act in a correlated 
manner. In this sense the system's behavior near its critical point indeed manifests in a 
macroscopic scale. 

Note traditional treatment of critical phenomena is classical in that the external parameters, 
temperature, magnetic field, etc, are classical. But it could also be quantum in nature, the 
system being driven by quantum noise or vacuum fluctuations. Quantum phase transition is 
brought about by the change of the coupling parameters when the system is at zero temperature. 
It results from the energy competition between different ground states as a function of the 
coupling parameter, and quantum entanglement in the system has been shown to provide a 
useful measure of such transitions. We hope to address this issue in a following paper when we 
turn our attention to the role of quantum entanglement in the description of MQP. 

Naively one would expect that at the critical point as the correlation length goes to infinity, 
the system would not depend on the microscopic details such as the lattice spacing in the Ising 
model, or the interaction amongst the microscopic constituents. But it does, and that makes it 
more interesting. It is this particular aspect of MQP which fascinates us, namely, what and how 
known microscopic details determine a quantum system's macroscopic attributes, or, even more 
challenging, how to decipher a system's microscopic features from the macroscopic phenomena. 
An example is the dialectic relation between emergent and quantum gravity pertaining to the 
micro and macro structures of spacetime [7] . 

On this point it is perhaps worth making the following observation: systems in a critical 
state can exist for a long time to certain classes of observers with specially interesting physical 
consequences, and many qualitative features in the ensuing phase depends on what came out 
of the transition. An example is inflationary cosmology where it is believed that the universe 
underwent a first (old) or second (new) order transition at the GUT (lO^^Ge^) scale, with a 



duration of 68 e-folding time to be able to account for the entropy content (10^*^ photons) in 
the present universe. During the exponential expansion the universe appears to be static for 
a local observer, its scale factor described by a de-Sitter -Einstein solution. In fact, in such 
an epoch between the metastable false vacuum and reheating to the true vacuum, one can use 
the language of critical phenomena to describe its dynamics, replacing the evolution by scaling 
(see, e.g., [8]). The microscopic features of the inflaton field - the quantum fluctuations which 
existed in the early universe - are magnified during this phase transition into galaxies of today. 
The effect of inflation is like a giant zoom lens at work [S] and the universe we observe today 
is in this sense a truly macroscopic quantum phenomenon [XU]. The exponential red-shifting of 
outgoing waves from a black hole which accounts for the thermal feature of Hawking radiation 
has a similar nature, acting like a microscope in a pictorial description [llj . 

1.3. Correlation in relation to fluctuations and coherence. Organization of this paper 
Unlike the situations of large A^ where an expansion in inverse N yields meaningful physical 
results, e.g., the leading order corresponds to a semiclassical limit, quantum correlations are less 
direct or transparent in their physical meanings. Correlation is related to fluctuations (as one 
form of fluctuation-dissipation relation indicates) and at the quantum level, quantum correla- 
tion can be a measure of quantum coherence (as the correlation history approach to quantum 
decoherence shows [12j. It is also related to but distinct from quantum entanglement (see work 
of Cirac [13]) It would be very useful to identify the conditions whereby these three properties 
could offer a measure (no matter how tenuous) of the macroscopic features of a quantum system 
and how they signify the quantum (albeit residual) nature of a macroscopic system. For this we 
use three examples to explore these aspects. Part I of this paper consisting of Sections 2 and 3 
discusses the relation between the loop, large A^ and nPI (n-particle irreducible) expansions, the 
first two signifying the appearance of quantum and macro features and how quantum correlations 
is placed in regard to them. In Sec. 2 we use the existence of a H theorem at the next-to-leading- 
order in the quantum mechanical 0{N) model to demarcate which order of large N expansion 
will entropy generation first ensue and thus can be (approximately) regarded as acquiring some 
macroscopic attributes. In Sec. 3 we provide a concrete example of the above connections in 
the Bose-Hubbard model for A^ atoms in an optical lattice. The nonequilibrium dynamics of 
this system is described by the equations of motion obtained from the closed-time-path (CTP) 
two particle irreducible (2-P/) effective action. We compare results from a large A^ expansion 
where M is the number of fields (see Eq. (j3.5p below) under second order coupling with the 
exact (numerical) solution to see how quantum correlation plays out against A/", the number 
of components of the N body wave function. Although M or N are not exactly indicators of 
the macroscopic, they offer a probe into the interplay between these three aspects: Correlation, 
fluctuations and coupling strength. Part II of this essay discusses the effective field theory (Sec. 
4) and effective infrared dimensional reduction (Sec. 5). Though we didn't dwell on this point 
as we believe it is quite well known, we want to point out the fact that the critical behavior of 
interacting quantum systems can provide a sense of macroscopia quite different from the sim- 
plistic depictions (referring to the size of the system or how many components it possesses). Key 
to this is the fact that coupling strengths can vary ('run') with energy scale according to the 
RG equations. In Sec. 4 we describe effective field theory, where this key notion well adopted 
in critical phenomena studies is put to good use in that the low energy phenomenology (or 
long wavelength behavior) is effectively independent of or largely insensitive to the high energy 
processes. In Sec. 5 we discuss the conditions where an interacting quantum system near the 
critical point may take on a lower-dimensional appearance. This effective infrared dimensional 
reduction is shown by an eigenvalue analysis: wherever the eigen-spectrum of the invariant op- 
erator of the order parameter field possesses a band structure with a gap separating the lowest 
mode or band from the higher sector, the infrared behavior of this system behaves effectively as 



that of a lower dimension. These examples serve to show the intricacy of even a simple notion 
such as 'macroscopic' for a quantum system, that it is invariably tethered with the interaction 
strength, enters into the correlation lengths and undergoes qualitative changes when the system 
is close to the condition of criticality. The characterization of quantum macroscopic phenomena 
by quantum entanglement is left for a later paper. 

Part I: Correlation hierarchy, nP/, A^LO-large N and H-theorem 

2. Correlation entropy and H theorem at NLO 

To prepare for our discussions on how quantum correlations play out in the quantum - macro 
issues we first introduce the nPI effective action with n=2 as example, which is a special case of 
the so-called master effective action [Ijj (MEA) , with the specific model of a quantum mechanical 
0{N) model (QMON) which we have used in our first paper (MQPl) to illustrate the large A^ 
expansion. We then recount (see [151 116] ) how nPI is related to the loop expansion on the 
one hand, where one customarily takes as a measure of quantum features, and the large A^ 
expansion on the other, where one may regard as signifying macroscopic features. (These latter 
two relations were expounded in MQPl.) We then discuss the implications of the existence of a 
H theorem proven in [T7] using the QMON model. In the next section we will use the example 
of the nonequilibrium dynamics of BEC atoms to illustrate how this threefold-relation amongst 
1/A^, 2PI and expansion in the (second order) coupling strength pertain to quantum and macro 
issues. 

We paraphrase the results from [17j where a correlation entropy for an interacting quantum 
field is proposed, and a proof of the existence of a H-theorem for the quantum mechanical 
0{N) model [181 US] is provided. For the former, Calzetta and Hu followed the paradigm of 
Boltzmann-BBGKY |2U| and proposed a correlation entropy (of the nth order) for an interacting 
quantum field [21] [22] obtained by 'slaving' the higher (n+1 th) order correlation functions in 
the Schwinger-Dyson system of equations (See [TSldB] and below). They then derived the closed 
time path (CTP) [2H] two particle irreducible {2PI) [23] effective action {EA) for this model 
up to the next-to-leading order {NLO) in 1/A^ [T^ and used the CTP 2PI EA to prove an 
H-theorem for the correlation entropy of a quantum mechanical 0{N) model at the NLO level. 
We begin with an introduction to the correlation hierarchy. 

2.1. Boltzmann-BBGKY Hierarchy and Schwinger-Dyson Equations 

As is well known in nonequilibrium statistical mechanics, truncation of the BBGKY hierarchy 
at a finite nth order yields a closed system [j. The system of equations describing the n-particle 
distribution functions are time reversal invariant. When a causal factorization condition such 
as Boltzmann's molecular chaos assumption is imposed - that is, assuming that the n+1 th 
correlation function can be factorized into a product of n th correlation functions initially but 
not finally (after the collisions) - time-irreversibility appears and an ff-theorem is obtained. 
This type of coarse graining of the hierarchy i.e., truncation plus causal factorization into a set 
of coupled equations for the n-particle distribution functions is called 'slaving' in the language 
of [15[I16]. Slaving of the n+1 -particle distribution function renders an otherwise closed system 
(of the n+1 th order correlation functions) an effectively open system (of nth order correlation 
functions) and ushers in the appearance of dissipative dynamics l22j. Noise and fluctuations 
iB] should also appear, as required by the fluctuation-dissipation relation (FDR) [27], now 



^ The truncation of the correlation hierarchy is not just an arbitrary mathematical procedure, it reflects the 
fact that in realistic conditions, measurement are of finite accuracy associated with limited resolution of the 
instruments. Thus the relevant physical degrees of freedom are often limited to the lower end of the hierarchy, 
namely, the mean field and two-point function. 



manifesting not just for an open system near equilibrium (as depicted in linear response theory), 
but for an effectively open system [151 IS]- 

A familiar example is the Boltzmann equation for dilute gases. At sufficiently low density the 
description of a molecular gas by a truncated BBGKY hierarchy with only the one-particle and 
two-particle distribution functions may be justified. Because of the low density one may further 
assume that before the collisions the colliding partners are uncorrelated, i.e., the two-particle 
distribution function factorizes into a product of single-particle distribution functions. This is 
the physical basis for Boltzmann's imposition of the molecular chaos hypothesis. Of course 
the colliding partners become correlated after interaction. This assumed distinction between 
initial and final conditions is the origin of the macroscopic arrow of time, the appearance of 
dissipation and foe entropy generation in the Boltzmann paradigm. FDR calls for a rightful place 
for fluctuations in the Boltzmann equation, as demonstrated in the derivation of a stochastic 
Boltzmann equation |28[ [29] . 

There are other approximation schemes that allow for the derivation of the kinetic equations 
from the correlation hierarchy [20j. To see the relation between correlation and interaction in 
this context of kinetic theory consider a system in which particles interact via a potential of the 
form V{ri,rj) = \v{ri,rj). For weakly coupled systems, i.e. / <^ 1, the hierarchy of equations 
can be written for a certain order of A. 

To 0(A) the equation for the one-particle distribution function, namely, the Vlasov equation, 
is closed, i.e., it does not involve any two-particle distribution function. We emphasize that in 
this case truncation alone not followed by slaving is sufficient to obtain a closed equation for 
the one-particle distribution function. There is no collision term for short-ranged interaction, 
but interactions between particles do exist represented by an average potential of all other 
particles present. Vlasov equation is useful for describing the kinetics of systems containing 
many particles interacting under long range forces. Recall we have seen it in our first paper 
as the equation obtained from the leading order large N expansion in the 0(N) model. Vlasov 
equation is obtained as a mean field theory with reversible dynamics and no H-theorem. 

To O(A^) the set of equations for the one-particle and two-particle distribution functions is 
also closed. By slaving the two-particle distribution function to the one-particle distribution 
function the Landau equation is obtained. Similar to the Boltzmann equation. Landau equation 
also has a collision term accounting for the collisions between individual particles at O(A^). 
Landau equation describes irreversible dynamics and its solutions satisfy the H-theorem. What 
is different from the Boltzmann equation is the form of the collision integral. This difference 
can be traced back to the absence of a term in the equation for the two-particle distribution 
function at O(A^) compared to the low density approximation. 

The situation in quantum kinetic field theory is completely analogous. By quantum kinetic 
field theory (see, e.g., [2T1 [22]), we are referring to the hierarchy of coupled equations for the 
relativistic Wigner function and its higher-correlation analogs, obtainable from the variation of 
the master effective action |15[ [TC] whose variation yields the Schwinger-Dyson equations. This 
is a quantum analogue of the BBGKY hierarchy, expressed in a representation convenient for 
distinguishing between microscopic (quantum field-theoretic) and macroscopic (transport and 
relaxation) phenomena |j. 

One may choose to work with a truncation of the hierarchy of the Wigner function and its 
higher correlation analogs, or one may instead perform a slaving of, for example, the Wigner- 
transformed four-point function, which leads (within the context of perturbation theory) directly 
to the relativistic Boltzmann equation [22] and the usual //-theorem. The truncation and 
subsequent slaving of the hierarchy within quantum kinetic field theory can be carried out at 

^ It should be pointed out that in order to identify the relativistic Wigner function with a distribution function for 
quasiparticles, one must show that the density matrix has decohered, and this is neither guaranteed nor required 
by the existence of a separation of macroscopic and microscopic time scales [30) . 



any desired order, as dictated by the initial conditions and relevant interactions. 

2. 2. Entropy from Slaving the Higher Correlations: The CTP 2PI EA 

A general procedure has been presented for obtaining coupled equations for correlation functions 
at any order / in the correlation hierarchy, which involves a truncation of the master effective 
action at a finite order in the loop expansion \lb\ 116] . By working with an I loop-order truncation 
of the master effective action, one obtains a closed, time-reversal invariant set of coupled 
equations for the first / + 1 correlation functions, 0, G, C3, . . . , Q+i. In general, the equation of 
motion for the highest order correlation function will be linear, and thus can be formally solved 
using Green's function methods. The existence of a unique solution depends on supplying 
causal boundary conditions. When the resulting solution for the highest correlation function 
is back-substituted into the evolution equations for the other lower-order correlation functions, 
the resulting dynamics is not time-reversal invariant, and generically dissipative. Thus, as was 
described before, with the slaving of the higher-order (Wigner-transformed) correlation function 
in quantum kinetic field theory, we have rendered a closed system (the truncated equations for 
correlation functions) into an effectively open system. In addition to dissipation, one expects 
that an effectively open system will manifest noise/fluctuations, as shown in [151116] for the case 
of the slaving of the four-point function to the two-point function in the symmetry-unbroken 
A^"^ field theory. Thus a framework exists for exploring irreversibility and fluctuations within 
the context of a unitary quantum fleld theory, using the truncation and slaving of the correlation 
hierarchy. The effectively open system framework is useful for precisely those situations, where 
a separation of macroscopic and microscopic time scales (which would permit an effective kinetic 
theory description) does not exist, such as is encountered in the thermalization issue. 

As a particular coarse graining measure the slaving of higher correlation functions to lower- 
order correlation functions within a particular truncation of the correlation hierarchy has 
several important benefits. It can be implemented in a truly nonperturbative fashion. This 
necessitates a nonperturbative resummation of daisy graphs, which can be incorporated in the 
truncation/slaving of the correlation hierarchy in a natural way. 

In [16J the authors are interested in the growth of entropy due to the coarse graining of the 
correlation hierarchy by slaving a higher correlation function. The simplest nonperturbative 
truncation of the Schwinger-Dyson equations for the A$^ field theory which contains the time- 
dependent Hartree-Fock approximation is the two-loop truncation of the master effective action, 
in which only the mean field (j), the two-point function G, the three-point function C3 are 
dynamical. All higher order correlation functions obey algebraic constraints, and can thus be 
expressed in terms of the three dynamical correlation functions. 

While this truncation of the Schwinger-Dyson equations is well-defined and could in principle 
be solved, it is disadvantageous because, as stated above, without some coarse graining, the 
system will not manifest irreversibility and will not equilibrate. Therefore we slave the three- 
point function to the mean field and two-point function, and thus arrive at an effectively open 
system. In principle, a systematic analysis of the coarse-grained dynamics of the mean field and 
two-point function should include stochasticity [151 [16]. 

The 2PI formalism is also suitable for addressing this question because, provided an auxiliary 
field is cleverly introduced, the 2PI CTP effective action can be found in closed form at each 
order in 1/N \il\ 131] . We now use the familiar quantum mechanical 0{N) model to illustrate 
how to derive the 2PI effective action and then explore the conditions for the existence of 
an H-theorem. This model was used in our first paper p] for the discussion of large N and 
macroscopia. 



2.3. 0{N) \X^ Theory 

As we recall the system dynamics is described by the Hamiltonian with variables Xa and their 

conjugate momenta P , where A,B are the 0{N) group indices, with 



F = ^ i^P^'p'' + M^XbXb + -^ {XbXbY 



(2.1) 



The classical action 



S 



dt ^ IxbXb - M^XbXb - -^ {XeXef 



(2.21 



where Mq and Aq are the mass parameter and coupling constant. We rescale Xb = v Nxb 



S = N dt -I xbxb - M'^xbxb - j (xb^b)^ 
Discarding a constant term, we may rewrite the classical action as 



S = N I dt - <^ XBXB 



M2 Vx 
Tx + ~"^"^ 



(2.3) 



(2.4) 



To set up the 1/A^ resummation scheme, it is customary to introduce the auxiliary field Xi 
writing 




XBXB 



M2 VA 



+ 



M 






(2.5) 



whence 



S = N I dt < -XBXB - X 



M^ XBXB 

~^ 2 



+ ^^ 



(2.6) 



From now on, we consider x ^-i^d xb as fundamental variables on equal footing. 

Because of the 0{N) symmetry, the symmetric point must be a solution of the equations of 
motion. For simplicity, we shall assume we are within this symmetric phase, and treat xb as a 
quantum fluctuation. We also split the auxiliary field X = X + X iiito a background field x and 
a fluctuation field x- The action becomes 



S — Sq + Si + S2 + s^ 
Sq is just the classical action evaluated at xb = 0, x = X'- 



So = jJdt i^x^-M^x 



(2.7) 



(2.8) 



Si contains terms linear in x and can be set to zero by a suitable choice of the background field 

X- 



5l = y|rft {X-M^}X 
S2 contains the quadratic terms and yields the tree - level inverse propagators, 

I. . X 1 



^2 



,~,2 



N I dt { -XBXB - -XBXB + r^X 



(2.9) 



(2.10) 



Finally 5*3 contains the bare vertex 
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To write the 2PI CTP EA we double the degrees of freedom, incorporating a branch label 
a = 1,2 (for simplicity, if not explicitly written, we assume that the label a includes the time 
branch, i.e., x = x {to))- We also introduce propagators G ' for the path ordered 
expectation values 



G 



Aa,Bb 



X^'^X^') 



(2.12) 



and F"^ for 



j,ab ^ /^a^b 



(2.13) 

Because of symmetry, it is not necessary to introduce a mixed propagator, for (x"'x^^) = 0. 
The 2PI CTP EA reads 



So 
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Sn 



X 



+ i J dudv \DAa,Bb{^^ V)G'^^'''\U, V) + yCabd{u, v)F^'' {u, v) 
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Tr InG + Tr \nF] + VQ 



(2.14) 



where, if the position variable is explicit, cn 



-C22 = 1, Ci2 = C21 = 0, 



DAa,Bb{u,v) = N6aB Cabdl - CabcX" ^{u,v), 



(2.15) 



and Cabc = 1 when all entries are 1, Cabc = ~1 when all entries are 2, and Cabc = 
otherwise. When we use the compressed notation, it is understood that Cab = CabS{ta,ti,) and 
Cabc = Cabc^{ta,tb)S{ta,tc)- Fq is the sum of all 2PI vacuum bubbles with cubic vertices from 
5*3 and propagators Q '^^ and F""^ . Observe that Tq is independent of x"^. 

Taking variations of the 2PI CTP EA and identifying x^ = x^ = Xi we find the equations 
of motion 



where Dab {u, v) = Cab 
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We shall seek a solution with the structure 
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ab 
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(2.16) 

(2.17) 
(2.18) 

(2.19) 
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P^B^-\u,v) 



(2.20) 



which is consistent with vanishing Noether charges. Then it is convenient to write 

pah^^jjab. Y[A^^Bb = 5ABPah; liab {x,y) = NQ^b {x,y) (2.21) 

The equations become 



Dab-'^ 

\cab - i 



G-1 



+ ^Pab = Q (2.22) 

ab I\ 



H-' 



ab 



+ Qab = (2.23) 



^[x{t)-M'}-^G'\t,t) = (2.24) 



Observe that 



_2«Q_ _ 2 STq 

These are the exact equations we must solve. The successive 1/A'^ approximations amount to 
different constitutive relations expressing Pab and Qab in terms of the propagators. 

For power counting in orders of N, observe that in any given Feynman graph each vertex 
contributes a power of N, each internal line a power of N^'^, and each trace over group indices 
another power of A^. We have both G and H internal lines, but the G lines only appear in closed 
loops. On each loop, the number of vertices equals the number of G lines, so there only remains 
one power of A'^ from the single trace over group labels. Therefore the overall power of the graph 
is the number of G loops minus the number of H lines. Now, since we only consider 2PI graphs, 
there is a minimum number of H lines for a given number of G loops. For example, if there are 
two G loops, they must be connected by no less than 3 H lines, and so this graph cannot be 
higher than NNLO. A graph with 3 G loops cannot have less than 5 H lines, and so on. 

We conclude that Tq vanishes at LO, and therefore Pab = Qab = 0. There is only one NLO 
graph, consisting of a single G loop and a single H line. This graph leads to 

r^^o = i-ih) (^-^) (- ^) ' 2iV (^ A^ ' CabcCdef I dudv H'^'' (u, v) G'^ {u, v) G'f (u, v) 

(2.26) 
Therefore, we get 

Pab = inCacdCbefH''G''f (2.27) 

Qab = irCacd.CbefG'^'^G •^ (2.28) 

The 2PI EA yields equations of motion for the (arbitrary time) two-point functions of the 
theory. Given a solution of these equations, in principle we may find the expectation values of 
a large family of composite operators at any given time. Suppose we adopt a coarse-grained 
description where we choose a certain number of these expectation values as the relevant variables 
to describe the system. Then there will be a single density matrix which has maximum von 
Neumann entropy with respect to the class of states reproducing the preferred expectation 
values. This maximum entropy density matrix is the reduced density matrix for the system, and 
its entropy its correlation entropy. The H theorem is the statement that the correlation entropy 
grows in time, when the correlations themselves are evolved using the equations derived from 
the 2PI EA truncated to some order in the 1/N expansion. For details in the proof refer to 



2.4-- H theorem at NLOLN and its implications 

It is clear that if we could solve the exact evolution, the H theorem would be manifest: Given 
/O (to) at the initial time to, solve the exact Liouville equation up to a time fi. Let p{ti) be 
the result. We extract the new expectation values from p{ti) and use them to construct the 
new maximum entropy density matrix p{ti). Then S[p{ti)\ > S[p{ti)\ , by definition, and 
•S" [p (ti)] = S [p (to)]) because the exact evolution is unitary, thus the H theorem. 

What we need to determine is whether, given the approximate dynamics for the expectation 
values provided by the 1/N scheme, the H theorem still holds. From our first paper (MQPl) 
in this series PP we know that to LO it does not because LO large N yields a mean field 
theory whose equation of motion, the Vlasov equation, is unitary. Thus there is no net entropy 
production and no //-theorem. In such a case, if one works with a Fock representation, for 
boson fields, the number of particles increases with time and can be used as a measure of field 
entropy, what Hu and Pavon [32] called an intrinsic measure of field entropy. (See also [33j.) 
However, this should not be confused with the correlation entropy of the Boltzmann kind under 
study for which the //-theorem is defined. 

Of course, we expect better approximations to be closer to the exact dynamics, therefore 
requiring less external control of (or rather, tampering with) the system. This reduces the 
entropy loss to the environment. Physically, whenever some higher order correlations are ignored, 
we expect to obtain an H theorem because such a description of the system is incomplete. This 
is the physical origin of correlation entropy, as explained in the Introduction. In order to have 
the different correlations evolving according to the 2PI 1/N equations appropriate to the desired 
order, instead of the accurate description from the full Schwinger - Dyson hierarchy, we must keep 
nudging it to conform to this artificially created condition due to our inability to comprehend the 
complete picture. This is the source of a new kind of noise arising from ignoring the higher order 
correlation functions, called correlation noise, which, when incorporated in the Schwinger-Dyson 
hierarchy, yields the (nth-correlation order) stochastic SD equations, as explained in [16] . 

2.5. 2PI EA, H-Theorem, quantum- classical and micro-Macro 

We can now ask the question: How does quantum correlations get registered in the 2PI EA 

and how does the existence of H-theorem at the NLO expansion pertain to the quantum-Macro 

issue? 

First, how do the quantum features show up in the 2PI effective action? Only equations of 
evolution for the first two moments can be obtained from the 2PI EA. It is known that [33| the 
equations of motion for the first two moments are identical for classical and quantum systems. 
Quantum corrections arise from the way the higher order moments are treated implicitly in the 
effective action approach. Because of this implicit dependence the 2PI EA approach is not 
easy to work with in addressing the differences of quantum and classical dynamics. One could 
consider going to higher orders, at least 3PI, and see if the third moment equations possess the 
terms absent in classical dynamics at a given order in large N . 

The micro-Macro (m-M) demarkation issue is as difficult to tackle, if not more, as the 
quantum-classical transition issue. Going from micro to macro requires coarse graining. QMON 
model is a "democratic" systemQ, hence the coarse graining should be "democratic" as well. The 
2PI approach provides us with a reduced description in terms of the two lowest order Green's 
functions, which we shall refer to as the relevant variables. Higher order Green's functions are 
referred to as irrelevant variables. If one keeps track only of the relevant variables information 
about irrelevant variables is lost, by choice. If these variables are coupled and one wishes to see 

* Democratic' is a figurative description of a system of identical particles and 'autocratic'refers to a system 
containing a distinct member from the rest. This corresponds to the Boltzmann versus Langevin paradigms 
describing an effectively open (by correlation order, which depends on the level of precision of observation) versus 
an open (by ab initio designation) system. See, |35) . 



how the relevant variables (our system) is affected by the irrelevant variables (its environment) 
one needs to take into account the back-action of the irrelevant variables. This back-action 
included reduced system (in terms of e.g., a reduced density matrix) would show dissipative 
dynamics and a H-theorem will obtain for the correlation entropy. 

The question is whether the H-theorem holds for the dynamics obtained at a certain order 
of the \/N expansion. There is no H-theorem at the leading order. This can be understood 
by realizing that the LO is equivalent to Gaussian approximation, and in this approximation 
correlations of order higher than two have no dynamics. As a result the mean field and the 
two point function contain all the available partial information because of the approximation 
imposed. The fact that H-theorem holds at NLO shows that the NLO approximation is capable 
of capturing the back-action from the irrelevant variables which engenders dissipation in the 
dynamics of the relevant variables. Very coarsely, in the large N perspective the NLO may be 
viewed as providing us a witness to the micro-macro transition. Note this does not mean that 
NLO large A'^ is a good approximation in general. Actually we know from Mihaila et. al. [18] 
that it fails at long times although the H-theorem holds at long times, albeit only shown for 
a vacuum initial state [T7|. The existence of an H-theorem, hence irreversibility, is a way of 
arguing that we have gone beyond the micro-domain. 

In the QMON model A^ plays a dual role. On the one hand A^ is the number of degrees of 
freedom. On the other hand N appears in the denominator of the interaction term in the system 
Hamiltonian 12 . 1 [ The large N limit is usually associated with the number of components A^ — )• oo 
only, however for the QMON model it is the weak coupling limit as well. If we would like to 
make an analogy to the classical BBGKY hierarchy of section [2711 large N approximation would 
have the same effect as the weak coupling approximation. The fact that there is no H-theorem 
at LO but there is an H-theorem at NLO is analogous to the fact that Vlasov equation doesn't 
lead to entropy production, whereas the Landau equation does so. 

Entropy increase is a signifier of irreversibility and irreversibility is usually associated with 
macroscopia. LO approximation gives reversible dynamics, while at NLO irreversibility emerges. 
This fact appears to be in contradiction with the association of irreversibility with macroscopia, 
since LO theory is successful in describing systems with more components A^ than NLO theory. 
The apparent paradox can be resolved by appreciating the dual role of A^ in the QMON model. 
The above example shows the interplay between system-size and interaction strength in the 
emergence of macroscopic qualities in a system. The system-size alone is not enough to qualify 
the system as macroscopic. 

The work of Mihaila et. al. |18| reveals a critical value of A^ below which the NLO 
approximation does poorly (Ac = 18.6 for the parameters they chose for the simulations). 
As was discussed in MQPl this is a way of quantifying when A^ is actually large enough for a 
certain order to yield good results. This is a stronger condition than the validity of H-theorem 
(which holds at NLO for any A"), since it requires the dynamics to follow the exact one, rather 
than just one condition of nondecreasing entropy. We may combine these two results into the 
following interpretation: The H-theorem is necessary for a satisfactory reduced description of 
a macroscopic system. It holds at NLO. NLO is a reasonable approximation for the system 
dynamics for N > Nc- Hence for N > Nc the reduced description is consistent. 

3. Model study: 2PI, NLOLN, loop expansions in Nonequilibrium BEC dynamics 

We begin with a description of the well-known Bose-Hubbard model and state our goal in 
this illustrated example. We introduce the 2PI generating functional to construct the 2PI 
effective action (EA) r[(j),G] and Green's functions. We perform a perturbative expansions 
on T2 and define the various approximation schemes. We then derive the equation of motion 
from the CTP 2PI EA and discuss the results under each approximation scheme, starting with 
the Hartree-Fock-Bogoliubov (HFB) approximation followed by second order expansions which 



includes the 1/N expansion and then the full second order expansion. The results from these 
approximations are compared with the exact numerical answers so as to exhibit the effects of 
quantum correlations, the coupling strength and how the results depend on increasing A^. 

The dynamics of an ultracold bosonic gas in an optical lattice can be described by a Bose- 
Hubbard model where the system parameters are controlled by laser light. For a one dimensional 
lattice the Hamiltonian is: 

i it 

where ^j and ^j are the annihilation and creation operators at the site i which obey the 
canonical commutation relations for bosons. Here, the parameter U denotes the strength of 
the on-site repulsion of two atoms on the site i; the parameter Si denotes the energy offset of 
each lattice site due to an additional slow varying external potential such as a magnetic trap 
and J denotes the hopping rate between adjacent sites. Because the next-to-nearest neighbor 
amplitudes are typically two orders of magnitude smaller, tunneling to them can be neglected. 
We denote the total number of atoms by A^ and the number of lattice sites by I. Only a one- 
dimensional homogeneous lattice with periodic boundary conditions is considered. 

A dimensionless parameter that is convenient to describe the different regimes of H is the 
coupling strength A = NU/IJ. Different from a homogeneous system without a lattice where 
at zero temperature the superfluid fraction is always unity, the presence of the lattice changes 
the superfluid properties and even at zero temperature, the superfluid fraction decreases with 
the lattice depth. For strong-coupling strengths [3B] I > lent-, 



2N. 



lcr^t-J^[2N + I+^{2N + Iy + I% (3.2) 

it is known that the ground state undergoes a quantum phase transition from a superfluid to 
a Mott insulator. In the weakly interacting regime, / <C 1, where tunneling overwhelms the 
repulsion, to a good approximation quantum fluctuations can be neglected and the properties 
of the system can be described by replacing the operator on the lattice site i by a classical c 
number. It can be said that most of the atoms are in the zero quasimomentum state. In the 
intermediate regime 1 <C ^ <C Icritf^ the interactions between the bosons can be very strong but 
the ground state is nevertheless a superfluid. For these interaction parameters a self-consistent 
HFB-Popov theory gives a good description of the system. However, different from the weak 
interacting regime where the depletion of the zero quasimomentum state is very small and has 
little effect on the superfluid properties, in this intermediate regime, depleted atoms spread 
over the central part of the band and reduces the superfluid fraction. As interactions are further 
increased the depleted population completely fills the band and cancels the superfluid properties. 
The system reaches the Mott insulator regime, where atoms are localized at each lattice site 
and the eigenstates of the system are almost Fock states with vanishing number fluctuations per 
lattice site. The dynamics in the intermediate regime is the focus of study in [37], where the 
superfluid properties are important but quantum fluctuations cannot be ignored. 

3.1. 2PI EA r[0,G] 

The flrst requirement for the study of nonequilibrium processes is a general initial-value 
formulation depicting the dynamics of interacting quantum fields. The CTP or Schwinger- 
Keldysh effective action formalism f23] serves this purpose. The second requirement is to 
describe the evolution of the correlation functions and the mean field on an equal footing. The 
two particle irreducible {2PI) formalism [23] where the correlation functions appear also as 
independent variables serves this purpose. By requiring the generalized (master) CTP effective 



action |14j to be stationary with respect to variations of the correlation functions an infinite 
set of coupled (Schwinger-Dyson) equations for the correlation functions is obtained which is a 
quantum analog of the BBGKY hierarchy. The 2PI effective action produces two such functions 
in this hierarchy. 

The classical action corresponding to the Bose-Hubbard Hamiltonian (|3.ip is given in terms 
of the complex fields <I>j and $1 by 

(3.3) 
where, as before, i denotes the lattice position, J is the hopping rate and U is the interaction 
strength. We limit the analysis to the case when no external potential is present and include 
only nearest neighbor hopping. To compactify our notation we introduce '^"(a = 1,2) defined 
by 

$, = $1, $* = cl>2. (3.4) 

In terms of these fields the classical action takes the form 

sm = jdtY: {^-ha,^imdt^\{t) + Jaab^'i^,m\{t) - ^{a,,^i{t)^\{t)f^ , (3.5) 

where N is the number of fields (two in this case), and summation over repeated field indices 
a,b = (1,2) is implied, hab and aab are matrices defined as 

hab = ^i ^ Q j ^'^^"(10) ^^'^^ 

After second quantization the fields $" are promoted to operators. We denote the expectation 
value of the field operator or mean field by (p^it) and the expectation value of the composite 
field by Gfht,t'). Physically, |9i>f(t)p is the condensate population and the composite fields 
determine the fluctuactions around the mean field. 



m) = i'^fit)), (3.7) 

nG'}^{t,t') = {Tc^fm'dt'))-{<^t{t)){^'l{t')). (3.8) 

The brackets denote taking the expectation value with respect to the density matrix and Tq 
denotes time ordering along a contour C in the complex plane. 

All correlation functions of the quantum theory can be obtained from the effective action 
r[i?!), G], the two particle irreducible generating functional for Green's functions parametrized 
by ipiit) and the composite field G1j{t,t'). To get an expression for the effective action we first 
define the functional Z[J,K] [2l] as 

Z[J,K] = e'/'^^l-^'^l (3.9) 



where we have introduced the following index lowering convention 

Xa = aabX\ (3.10) 



The functional integral ()3.9p is a sum over classical histories of the field <1>^ in the presence of 
the local source Jja and the non local source ^ijab- The coherent state measure is included in 
D^. The addition of the two-particle source term is what characterizes the 2PI formalism. 
We define r[0, G] as the double Legendre transform of M^[J, K] such that 

OJiayt) 



Expressing J and K in terms of (j) and G yields 



r[0,G] = Ty[J,K]-y'(it^Ji,(t)0^(t)-^^dfdt'^</.,"(i)<^,^(t')K,,afe(t,t') (3.13) 

i ij 

-\ j dtdt' J2 Gijit^ t')Kijab{t, t'). 



From this equation the following identity can be derived: 



<5<A?(t) 



J^a{t) + J dt'J2{Kijad{t, t'))cp%t') j , (3.14) 



5T[(P,G] h 
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-K,jab{t,t'). (3.15) 



ij 



In order to get an expression for r[(p,G] notice that by using ()3.9p for VF[J,K] and placing 
it in ()3.13p for T[(j), G], it can be written as 

exp(^^r[,^,G]) = n/^^"exp|^(^5[$] + y'dt,J,,(t)[cI>^(t)-</>f(t)] (3.16) 

+1 J dUdt'j ($^(t)K,^,,(t, t')^]{t') - <Pt{m^Jab{t. t')<t^]{t')) - ^ Tr Gk) I 



^ J di,d£^ mt) - <i>m j^^^ [m') - 4{t') 



+ TrG- 



5G 



where Tr means taking the trace. For simplicity we have denoted J dtJ2i by / dti. Defining the 
fluctuation field, ipf = $f — (pf, we have 

T[<p,G] - TrG^-^^^ = -mnl[jD^-e^p(^^S[4>,G;^]^ (3.17) 

5[0,G;^] = S[^ + ^]-J dU ^-^^^tit) -\j dUdt', ^m /^^(f j) ^^it')- (3-18) 



By introducing the classical inverse propagator iD ^{(p) given by 



r\ -1 SS[(f> 



iA,a.(i,t')-' = ..„;":i..,. (3.19) 



{6ijhabdt + J{6i+ij + Si-ij)aab) S{t - t') 

U_ 



—^ {2(t)ia{t)(kib{t) + aabCpt{t)cpic{t)) Sij6{t - t'), 



the solution of the functional integro-differential equation ()3.17p can be expressed as 

T[^,G] = S[(l)] + ^TilnG-^ + ^Tr D-\(I))G + T2[(p,G] + const. (3.20) 

The quantity r2[0, G] is conveniently described in terms of the diagrams generated by the 
interaction terms in S[(f) + (/?] which are of cubic and higher orders in (p. 

Sint[(l^ + ^] = -^Jdti {^ib{t)v1{t)f -^Jdti ^tit)^ia{t)v'{t)ipib{t). (3.21) 

It consists of all two-particle irreducible vacuum graphs (the diagrams representing these 
interactions do not become disconnected by cutting two propagator lines) in the theory with 
propagators set equal to G and vertices determined by the interaction terms in £'[(/> + 93] . 

Since physical processes correspond to vanishing sources J and K, the dynamical equations 
of motion for the mean field and the propagators are found by using the expression (j3.20p in 
equations p.l4p and (j3.15p . and setting the right hand side equal to zero. This procedure leads 
to the following equations: 

habhdt<Ptit) = -J{^^+la{t) + (pi-lait)) + ■^(0^d(t)(/)f (t) + Gfi ,{t,t))<j)^a{t) + (3.22) 
^-.{{Guad{t,t) + G,,da{t,t))<Pf{t) '^r2[</',G'] 



^yy^naay^,^; , ^,,aa,., .,,,., ,., ^^^^^^ 



and 



^ija6(*'*) — Dijab{t,t) —T,ijab{t,t), (3.23) 

^sG-Ht,t'y 



i,ab{t,t') = 2i ^ J];^ J^ . (3.24) 



Equation (|3.23p can be rewritten as a partial differential equation suitable for initial value 
problems by convolution with G. This differential equation reads explicitly 

KndtG,f{t,t') = -J{^Gt^^{t,t') + Gti^{t,t')) + ^{<I^UM(t))G'^{t,t')+ (3.25) 

2JJ f 

j^^fit)Gt^it,t')Mt) + iJ d4j:utX)Gt,{t'\t') + id'^%,^^^ 

The evolution of </>" and G is determined by Eqs. (|3.22p and (|3.25p once r2[4>, G] is specified. 
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Figure 1. Two-loop (upper row) and three-loop diagrams (lower row) contributing to the 
effective action. Explicitly, the diagram a is what we call the double-bubble , b the setting-sun 
and c the basket-ball. 

3.2. Perturbative Expansion ofT2{(j),G) and Approximation Schemes 

The diagrammatic expansion of r2 is illustrated in the accompanying Figure where two and 
three-loop vacuum diagrams are shown. 

The dots where four lines meet represent interaction vertices. The expression corresponding 
to each vacuum diagram should be multiplied by a factor (— i)'(i)*~^ where / is the number of 
solid lines and s the number of loops the diagram contains. 

The action T including the full diagrammatic series for r2 gives the full dynamics. It is of 
course not feasible to obtain an exact expression for r2 in a closed form. Various approximations 
for the full 2PI effective action can be obtained by truncating the diagrammatic expansion 
for r2. Which approximation is most appropriate depends on the physical problem under 
consideration: 



3.3. 
(i) 



(ii) 



The standard approaches 

Bogoliubov (One-loop) Approximation: 

The simplest approximation consists of discarding r2 altogether. This yields the so called 
Bogoliubov or one-loop approximation whose limitations have been extensively documented 
in the literature ([38 } l39 ] ). 

Time-dependent Hartree-Fock-Bogoliubov (HFB) Approximation: 

A truncation of r2 retaining only the first order diagram in U, i.e., keeping only the double- 
bubble diagram. Fig. a, yields equations of motion of and G which correspond to the time 
dependent Hartree-Fock-Bogoliubov (HFB) approximation. This approximation violates 
Goldstone's theorem, but conserves energy and particle number |40^ I^Tl I42j ) . The HFB 
equations can also be obtained by using cumulant expansions up to the second order [43] 
in which all cumulants containing three or four field operators are neglected. The HFB 
approximation neglects multiple scattering. It can be interpreted as an expansion in terms 
of Ut/J, (where t is the time of evolution) and is good for the description of short time 
dynamics or weak interaction strengths. 



(i) Second Order expansion: 



A truncation retaining diagrams of second order in U containing besides the double-bubble, 
also the setting-sun and the basket-ball. By including the setting-sun and the Basket-ball 
in the approximations we are taking into account two particle scattering processes [44J. 
Second order terms lead to integro- differential equations which depend on the time history 
of the system. 

(ii) Large- AA approximation 

The l/M expansion is a controlled non-perturbative approximation scheme which can 
be used to study non-equilibrium quantum field dynamics in the regime of strong 
interactions |31j. In the large M approach the field is modeled by M fields and the quantum 
field generating functional is expanded in powers of l/M. In this sense the method is a 
controlled expansion in a small parameter but unlike perturbation theory in the coupling 
constant, which corresponds to an expansion around the vacuum, the large M expansion 
corresponds to an expansion of the theory about a strong quasiclassical field. 

In [37j numerical solutions to the equations of motion for a moderate number of atoms and 
wells up to the second order in the coupling constant U were given. These exact many body 
solutions are useful for determining the range of validity of the three types of approximations 
described above, namely, a) the HFB, b) the full second order and c) NLO large J\f expansion 
up to second order in U (the shorthands HFB, 2nd and 1/J\f respectively are used in the figure). 
For the effect of the total number of atoms on the dynamics Rey et al have given the numerical 
solutions for a double well system with fix ratio UN/ J = 4 and three different values of N: 
N=20, 40 and 80. The results are given for the evolution of the atomic population per well 
(Fig. 5), the condensate population per well and total condensate population (Fig. 6) and the 
quasi-momentum intensities ( Fig. 7). (To make the comparisons easier the numerical results 
obtained for the three different values of A^ are scaled by dividing them by the total number of 
atoms. In this way for all the cases we start with an atomic population of magnitude one in the 
initial populated well.) 

In the exact dynamics we see that as the number of atoms is increased the damping effects 
occur at slower rates. This feature can be noticed in the quantum dynamics of all of the 
observables depicted in the Figures 5 to 7 of [37J. The decrease of the damping rates as the 
number of atoms is increased is not surprising because changing the number of atoms alters the 
quantum coherence properties of the system. As shown in reference [45j, the collapse time of the 
condensate population is approximately given by tcoii ~ ^^-^ where a is the variance of the initial 
atomic distribution and trev is the revival time which depends on the detailed spectrum for the 
Hamiltonian. In the kinetic energy dominated regime t^ev ~ hN/J (see |46J), thus tcoii ~ — • For 
our initial conditions the variance is proportional to a = ^/N so tcoii ~ VN. Besides damping 
rates, the qualitative behavior of the exact quantum dynamics is not affected very much as the 
number of atoms is increased for a fixed UN/ J. 

The improvement of the 2PI approximations as A^ is increased, as a result of the increase in 
the initial number of coherent atoms is in fact observed in the plots. Even though the problem of 
underdamping in the HFB approximation and the overdamping in the second order approaches 
are not cured, as the number of atoms is increased, we do observe a better matching with the 
full quantal solution. The l/M expansion shows the fastest convergence, as can be more easily 
observed in the quasi-momentum distribution plots. Fig. 7. The better agreement of the l/M 
expansion relies on the fact that even though the number of fields is only two in our calculations 
the l/M expansion is an expansion about a strong quasiclassical field configuration, a feature 
we have explained in MQPl. 

Part II: Effective Scale and Infrared Dimension 



4. Effective field theory, coupling and the effective scale of a quantum system 

Effective field theory is a very powerful tool in quantum field theory, and in particular it gives 
a new point of view about the meaning of renormalization. 

Consider a quantum field theory with a characteristic energy scale £"0, and suppose that we 
are interested in physics at some lower energy scale E. The basic idea of effective field theory 
is that one can choose a cutoff energy scale A at or slightly below Eq, and divide the fields in 
the path integral into low- and high- frequency parts, (p = 4>l + fpH, where (pLi't'H) are the field 
modes with frequencies lower(higher) than A. Now integrate over the high frequency fields (pn 

D(l)HDct)Le'^^^"'^^'> = f D(l)Le''^^^'^^'> (4.1) 



where e^^^^'^'^' = ^ Dcjiue^^^'^^''^^' . The S\{cI)l) is called the low-energy or Wilsonian effective 
action. It's different from the IPX effective action which are generated from integrating over all 
frequencies but keeping only IPI graphs. 

One can expand the low-energy effective action in terms of local operators which are consistent 
with the symmetries of the problem, 

Sk{H) = j d^'xY^g^Qi. (4.2) 

i 

In units oifi = c = 1, if an operator Oi has unit E '■ where di is the dimension of operator Oi 
then the parameter gi has unit E ^^ where D is the spacetime dimension. In the case of a free 
scalar field in D dimensions, with action 



S= I d'^xi -d^4>d^<P - -m202 ) . (4.3) 



The action S is dimensionless, so the dimension of (/>(x) is determined from the kinetic term to 
be [(f)] = -^— and the dimension of m^ is [m?] = 2. An operator Oj constructed from M cp's and 
N derivatives has dimension di = M (— ^) -|- N. Define dimensionless couplings A^ = A'^'-~^gi. 
Now for a process at energy scale E, one can dimensionally estimate the contribution from the 
i-th term of the effective action 

d'^xgA ~ Hjf'-''. (4.4) 

If di > D, this term will become less and less important at lower and lower energies, and is 
called an irrelevant operator. If di < D, the operator is more and more important at lower and 
lower energies and is called a relevant operator. If di = D the operator is equally important 
at all scales and called marginal. In most cases there is only a finite number of relevant and 
marginal operators so the low energy physics depends only on a finite number of parameters. 

The low energy physics depends on the short distance theory only through the relevant and 
marginal couplings, and possibly through some leading irrelevant couplings if one measures small 
enough effects. 

It should be pointed out that the coupling strength measured will vary with energy (scale). 
This can easily be seen from the point of view of effective field theory. We could write the path 
integral as 

DcPke"^'-^^^^ = I Dc^fc/e^-^A'^"^*'). (4.5) 

fc<A Jk'<K' 

The action S\{(j)k) contains all momentum modes up to some maximum value A, whereas the 
new effective action S'j<^,{(/)k') contains momentum modes up to A' < A. If one take A' to be 



infinitesimally smaller than A, S'^, will be infinitesimally different from S\. This generates a 
differential equation, the Wilson equation 

^ = F{Sa), (4.6) 

where F is some functional of the action. The Wilson equation is the renormalization group 
flow equation in an infinite dimensional space and it gives the change in action as a function of 
cutoff. Since 

Sk{<P) = j d^xY^giQi, (4.7) 

i 

the Wilson equation then gives 

^ = n{a^})- (4.8) 

If one linearizes a solution of the equation, irrelevant operators will correspond to those directions 
with negative eigenvalues whereas the relevant operators corresponds to those directions with 
positive eigenvalues. 

5. Infrared behavior of quantum fields: Finite-size / Topology effects 

In the above we pointed out that the infrared or large scale behavior of a system depends on many 
factors: a) Interaction strength: difference between a free theory (Gaussian) and an interacting 
theory, and for an interacting theory the Gaussian fixed point as different from non-Gaussian 
ones, relevant versus irrelevant fixed points; b) Coupling strength which runs with energy via 
the RG equation. All this show the intricacy of the notion of 'large scale' which we need to bear 
in mind in describing 'macroscopic' quantum behavior. In this section, we demonstrate how 
the infrared behavior of a quantum system can be altered qualitatively by the geometry and 
topology of space, or more precisely, the 'finite size' effect. This can be seen from analysis based 
on rather general arguments, e.g., on whether there exist a band structure in the spectrum of 
eigenvalues of the invariant operator of the order parameter field. If so there is a dimensional 
reduction effect where the IR behavior of the system becomes that of a lower dimension. What 
this conveys is that, for a quantum system near the critical point where its correlation length 
increases to (near) infinity its large scale or IR behavior or macroscopic features can become 
sensitive to the size and topology of the underlying spacetime. Cosmological phase transition is 
an extreme yet valid, even realistic, manifestation of this feature. A more rigorous and complete 
treatment of this effect is given in Hu and O'Connor [17] from which our present exposition is 
adapted. 

Consider an A^-component self-interacting scalar field <l>"(a = 1, • • •,A^) on a manifold of 
dimension D, coupled to the background spacetime with curvature R and coupling constant ^ 
(conformal coupling for ^ = and minimal coupling for ^ = 1) described by the action 






(5.1) 



where M| = m^ + (1 — 0-^/6 and V^ = — ^fgd ^i_{^g^^ dp) is the Laplace-Beltrami operator for 
scalar fields. For convenience we develop the formalism here in Euclidean space version with 
signature (-|-,-|-,-|-,+). We begin with a decomposition of the field ^"^ into a background field cfP" 
and a fluctuation field (/?'* , i.e., $" = c^** -|- (/j". The background field 4'°' is required to satisfy 
the classical equations of motion with an arbitrary external source. Such a shift eliminates the 



linear term in the fluctuation field, which is equivalent to performing a Legendre transform. The 
resultant action is 



Mi*"" 



2 , T\j2 , ^I2^ xab , ^ Jajb 
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V^ + M^ + -(/)^ jS"" + - 



J) , -^ 7a, a,.2 , ^ ,4 



Z+g^VV + ^V^n- (5.2) 



In the Feynman diagrammatics there are two kinds of vertices: a four point vertex proportional 
to A and a three point vertex proportional to X(j)a{x). 

The effective action r[(^] is obtained by functionally integrating over the fluctuation fields: 

e-rM=y'[d(^]e-^[^'^]. (5.3) 

The wave operator A for the fluctuating field is given by 

A'^' = (V^ + Mf)^ + (V^ + Ml) U' - ^\ (5.4) 

where 

Ml = Ml + ^4>\ Ml = Ml + ^^K (5.5) 

Here (ffcf) /(j)^ and 5 —<j)"'(f) /(j)^ are orthogonal projectors, the former projects along the direction 
in the internal space picked out by 0" and the latter projects into an (A^-l)-dimensional subspace 
orthogonal to the direction of cf)"". Note that the operator V^ does not commute with the 
projectors unless (^" is a constant. 

When the direction in group space picked out by (j)"" does not vary from point to point around 
the manifold (this does not necessarily imply (j)^ is a constant) the Green's function for A is 
given by 

G.' = ft^ + G. (.--«!), (5.6) 

where the Gi (i = 1, 2) are the Green's functions for the operators V^ + Mf 

The one-loop correction is related to the sum of the logarithms of the determinants of the 
fluctuation operators. So in this case, by using the projection operators, the one-loop effective 
action is given by 

r[<^] = 5[0] - ^ Tr In Gi - ^^ Tr In G2 



1 ,^ _ ^^ N -l.^ 
2 



S[(l>] + -^di\n\n + ^—^di\n\2i, (5.7) 

where S\^ is the classical action, Aj/ and di are the eigenvalues and degeneracies of the operators 
V^ + Mf. 

The determinant of A is formally divergent and needs to be regularized. There are a number of 
commonly used regularization methods. If the space is Riemannian and has sufficient symmetry 
so that the spectrum of the invariant operator is known explicitly, then the ("-function method 
is probably most convenient. The generalized C function is defined by 

C(^) = Y.^I^-^^^)-\ (5.8) 

n 

where An are the eigenvalues of the operator A on the Euclideanized metric obtained by a 
Wick rotation to imaginary time t = it. Here a constant mass scale fi is introduced to make 



the measure d[ip\ of the functional integral dimensionless. Using the regularization method of 
Dowker and Critchley, one can express the one-loop effective potential as 

y(i)(0) = An{Vol)-\C'{Q)+m/v]. (5.9) 

The ultraviolet divergence in V^^' can be canceled by the addition of counterterms which is 
not our concern here. The regularized effective action is useful for the analysis of the infrared 
behavior of quantum fields in a curved spacetime or spaces of nontrivial topology or finite 
extent, as shown below. Note, however that one-loop results are insufficient. Leading infrared 
contributions in higher loop terms need be included by using composite operator techniques. 

We learn from the work of [U] that in spacetimes with some compact dimensions the lowest 
mode of the fluctuation operator has the strongest effect on the symmetry behavior of the 
system. When the lowest mode is massless it will give the dominant contribution to the effective 
action. The low energy behavior corresponds to a lower-dimensional system. In what follows 
we first give a formal derivation of infrared dimensional reduction by examining the result of 
the decoupling of the higher modes (or bands) in the functional integral for the effective action. 
We will then give a physical explanation in terms of the correlation lengths and the notion 
of effective IR dimension (EIRD). An alternative way of seeing this problem of dimensional 
reduction is by spectral analysis. This is applied to direct product spaces with some compact 
dimensions and to spaces which can be reduced to product spaces. 

5. 1 . Decoupling of the Higher Modes (or Bands) 

Let us examine systems where the eigenvalues of the fluctuation operator takes on a band 
structure. By band structure we mean that the eigenvalues occur in continua with each 
continuum having a higher lowest eigenvalue than the previous one. This is true for fields 
on spacetimes with compact sections or for fields with discrete spectrum (e.g. the harmonic 
oscillator) . The procedure is to expand the fields in terms of the band eigenfunctions and convert 
the functional integral over the fields to an integral over the amplitudes of the individual modes. 
On a manifold with topology R'^ x B^ where B is compact, consider quantum fields where 
the fluctuation operator A in (j5.4p has the general form of a direct sum of operators D and B 

A''\x, y) = D^^x) + B^''{y) (5.10) 

with coordinates x on R"^ and y on B^. Assume that the eigenvalues a;„ associated with the 
eigenfunctions i/jniv) oi B""^ are discrete: 

i?«''Vn(y)=t^f^n(y) (5.11) 

Decomposing the field ip"'{x,y) in terms of ipniu) 

^-{x,y) = J2^-^{x)My) (5.12) 

n 

one obtains for the quadratic part of the action 

1 J d'^xd^^^A'^V = \j d'^X^^J^mD'^'^l + C^^jVnn^V'^™] (5-13) 

where fnm = J d!'yipniy)ipmiy)- When ipn are properly normahzed /„„ = 6mn (we will make 
such a choice here) the resulting theory in terms of the new fields 995^ will involve massive fields 
with masses determined by the eigenvalue matrix A"^ of the operators V^ + Mj [see ()5.7p ]. even 



if the fields in terms of the old variables appeared massless. We will take the smallest eigenvalue 
to be given by n = and assume that its only degeneracy is labeled by the indices a and b. 
Assume also that the operator D is simply minus the Laplacian V^ on R times 6 , the n = 
mode is then governed by the action whose quadratic term is 

yd'x[^-i-Vl)6^'4 + u^'foo^^,4] (5.14) 

Thus this appears like a d-dimensional field with an apparent mass matrix ujq . For the case of 
an N component A0^ theory the action after this decomposition takes the form 

+ ■^allm'A^m + -^^hnim'/k'Pn'P\^m]^ (5-15) 

where 5-^^^ = ^ d!'yCp'''4)n'>Pi^m and fknim = ^ <f'y'>Pk'^nipi>'>Pm- The effective action is now given 
by the functional integral 

e-rM=|[d<^a]e-S[0+^]. (5.I6) 

The interesting case occurs when the lowest eigenvalue approaches zero. At low energy the 
Appelquist-Carazzone decoupling theorem assures us that with higher modes decoupled from 
the dynamics, the infrared behavior is governed by the lowest band. We are then left with a 
purely lower dimensional theory. The higher modes do play a role in the ultraviolet divergences 
present in the theory (e.g. renormalization problem in Kaluza-Klein theories) and therefore 
determine the high energy running of coupling constants, but the infrared region of the theory 
is governed by the lower d-dimensional theory. 

5.2. Correlation Length and Effective Infrared Dimension 

The above result of dimensional reduction from a formal derivation of mode decoupling can be 
understood in a more physical way by using the concept of effective infrared dimensions (EIRD). 
By effective IR dimension we mean the dimension of space or spacetime wherein the system at 
low energy effectively behaves. One well-known example is the Kaluza-Klein theory of unifica- 
tion and cosmology. There, starting with an 11-dimensional spacetime with full diffeomorphism 
symmetry, after spontaneous compactification it reduces at energy below the Planck scale to 
the physical four dimensional space with GL(4,R) covariance and a seven-dimensional internal 
space with symmetry group containing the standard SU-^ x SU2 x Ui subgroups of strong and 
electroweak interactions. For observers today of very low energy the effective IR dimension of 
spacetime is four, even though the complete theory is eleven dimensional. By the same token, 
Einstein's theory of general relativity can presumably be regarded as the effective IR limit of an 
otherwise more complete and fundamental theory of quantum gravity, such as the superstring 
theory. For curved-space symmetry breaking considerations, the EIRD which the system "feels" 
is governed by a parameter r/ which is the ratio of the correlation length H and the scale length L 
of the background space rj = H/L. For compact spaces like S"^, L is simply 27r times the "radius" 
of 5^, the only geometric scale present. For product spaces R'^ x B^ with some compact space 
B, there are two scale lengths: L^ is finite in the compact dimensions and L,^ = 00 in the non- 
compact dimensions. Examples are Kaluza-Klein theories d=4, b=l, 6, 7, finite temperature 
field theory (imaginary time formalism, L = /3 inverse temperature) d=3, b=l, and Einstein or 
the spatially-closed Robertson- Walker universe d=l, b=3. 



The symmetry behavior of the system (described here by a Xip'^ scalar field as example) is 
determined by the correlation length H defined as the inverse of the effective mass M^fj related 
to the effective potential Veff by (we use subscript eff to denote quantities including higher loop 
corrections) 

n O Vpf f , o ^curvature— induced mass M? „ . , . 

^^ ^J,2 '4>min fi// \+ radiative corrections ' ^ ' 

It measures the curvature of the effective potential at a minimum energy state ((^ = for the 
symmetric state, or the false vacuum, (p = (j)min for the broken-symmetry state or the true 
vacuum.) The effective mass is defined to include radiative corrections to the same order cor- 
responding to the effective potential. (This quantity is called the generalized susceptibility 
function in condensed matter physics.) The critical point of a system is reached when H — t- c« 
or Mgff — )• 0. In flat or open spaces or for bulk systems, the critical point can be reached with- 
out restriction from the geometry (note that in dynamical situations, exponential expansion 
can effectively introduce a finite size effect equivalent to event horizons, see |l8j). However, 
in spaces with compact dimensions, the correlation length of fluctuations can only extend to 
infinity in the remaining non-compact dimensions, and thus the critical behavior becomes ef- 
fectively equivalent to a lower d-dimensional system. One can also think of H as the Compton 
wavelength A = In/Meff of a system of quasi-particles with effective mass M^fj- Any fine 
structure of the background spacetime with scale L is relevant only if A < L. Thus when A is 
small or ?7 <C 1, (far away from critical point, at higher energy, with higher mode contributions) 
it sees the details of a spacetime of full dimensionality. At this wavelength, the apparent size 
of the universe is large in both compact and non-compact dimensions. When A — )• oo or 77 S> 1 
(near critical point, IR limit, lowest mode dominant) finer features in the compact dimensions 
will not be important. The apparent 'size' of the universe will be measured by the non-compact 
directions and the EIRD is measured by the number of non-compact dimensions. The value of 
T] getting very large is an indication of when dimensional reduction can take place. 

Notice that in flat space {R'^) critical phenomena the effective potential V^ff (free energy 
density) depends on the coupling constants of fields which run with energy and temperature. 
In curved-space coupling parameters run also with curvature or the scale length of the space. 
This makes the concept of EIRD even more interesting, as there is now an interplay between H 
and L, and rj can either decrease or increase with curvature. For example, for A</)^ fields in the 
Einstein universe near the symmetric state = 0, the EIRD is equal to 1 but near the global 
minimum of broken symmetry state cf) = (pmin is equal to 4. Near the symmetric state, r] ^ 1 
signifies reduction of EIRD to one. This is consistent with the theorem of Hohenberg, Mermin 
and Wagner (for statistical mechanics on a lattice) and Coleman (for continuum field theory) 
which states that in dimensions less than or equal to two, the infrared divergence of the scalar 
field is so severe that there could be no possibility of spontaneous symmetry breaking: the only 
vacuum expectation value for <j) allowed is zero. Away from the region ~ the one-dimensional 
behavior no longer prevails. Indeed a global minimum of the effective potential exists at (pmin- 
Near (pmin, ?? ^ 1 and decreases with curvature. Thus the apparent size of the universe near the 
global minimum actually increases with increasing curvature. There is therefore no dimensional 
reduction and the system has a full 4-dimensional IR behavior. A transition to the asymmetric 
ground state is not precluded as symmetry breaking via tunneling is in principle possible. The 
complete picture extending from (f) = to (p = (pmin is a combination of one dimensional and 
four dimensional infrared behavior. Similar arguments can be applied to other spacetimes or 
field theories. Using this notion one can understand, for example, why it is often said that at 
high temperatures (small radius limit of 5"^) the finite temperature theory becomes an effective 
three-dimensional theory. 



5.3. Dimensional Reduction: An Eigenvalue Analysis 

In the above we have introduced the notion of effective infrared -dimension and suggested the 
ratio of the correlation length H to the geometric scale of spacetime L as a measure of the 
conditions for the system to behave effectively as in a reduced dimension in the infrared regime. 
We suggested that for product spaces R x B with some noncompact dimension d, the effective 
IR dimension is usually just d. We will now verify this assertion by analyzing the spectrum of 
the fluctuation operators in these spacetimes directly. 

Consider for simplicity direct product spaces. Examples are a) physical cosmological space- 
times with topology R^{time) x S^{or R^,H^,T^), b) Kaluza-Klein cosmology with M^ 
(Minkowski) xS^ (or other internal space), c) finite-temperature (imaginary-time) theory has 
R'^ X S^. Take a simple example S'^ x S^ for illustration. (This could be the spatial geometry 
of a "handled" Gowdy universe.) Similar reasoning can be extended to a wide range of product 
spaces. 

From (j5.7p the wave operator A governing the fluctuation fields in the large A^ limit is 
A = V^ + M2, where V^ is the Laplace-Beltrami operator on a general curved spacetime. For 
S'^ X 5^ with radii 02 and oi respectively, V^ is a sum of the total angular momentum operator 
L on S^, and L^ on S^ 

V2 = :^ + M (5.18) 

The eigenfunction is a product of Y£m{0, (p)e^^^ belonging to the eigenvalues 



a 



where n = {i,m,n),i = • -oOjm = —i to i,n = it 1,- • • The eigenvalues of A are then 
An = «^n + -^2 • ^^ ^^^ infrared region, we are interested in the contribution of the lowest eigen- 
value (zero mode) to the effective potential. We will consider the two limiting cases of 1) S'^ x R^ 
and 2) R^ x S^ obtained when ai and 02 — )■ 00 respectively, and show that the EIRD is equal 
to 1 and 2 respectively. 

The effective potential 14// can be constructed from the (^ function ({0) and its derivative 

cm by dssD. 

C(-) = E(^"'^n)-^ = ^^^ E(2^ + l)[^^ + ^ + Mir (5.20) 

where 2i + l is the degeneracy of m. In the limiting cases the summation over discrete quantum 
numbers will be replaced by integrals of the form / k dk. We will derive the dimensionality 
of the reduced system by finding D. Thus in case a) ai — )• c«, the lowest eigenvalues belong to 
the band £ = Oasai— t-oo, k = n/ai assumes continuous value and (^{i') becomes 

/oo 
dkik"^ + Ml)-", ai -^ 00, M2a2 << 1 (5.21) 

-00 

This is a one-dimensional integral. In case (b) 02 — )• 00, the lowest eigenvalues come from the 
lowest band n=0 given by the first term in (j5.19p . where i assumes continuous values as 02 — )• cxd. 
Now call k = i/a2,C,{^) becomes (n = 0) 

/•oo 

C{v) ~ 2a2 / dkkik^ + M'^y , 02 -> 00 , M2ai « 1 (5.22) 

Jo 



The extra factor of k comes from the degeneracy of m for nonzero i. The integral is thus 
2-dimensional, as expected. 

Physicahy, the two hmiting cases may represent two different symmetry states of the system: 
one with cyhndrical symmetry where the quantum states with good quantum number m are 
eigenstates of the L^ operator, the other with spherical symmetry with good quantum numbers 
i associated with the L^ operator (on S'^). As a result of symmetry breaking the original system 
can end up in one of these states with a different symmetry. External perturbations with a 
particular symmetry and sufficient magnitude will influence the selection of the end state of 
the system. A simple example in elementary quantum mechanics is the Zeeman versus the 
Paschen-Bach effects associated with many electron atoms subjected to an external magnetic 
field. The total magnetic quantum number will become a good quantum number at very strong 
fields. More sophisticated examples can be found in the consideration of classes of internal spaces 
admitting larger symmetry groups with the strong and electroweak gauge groups as subgroup 
in the Kaluza-Klein theory. By the same analysis, it is not difficult to see that the EIRD of any 
homogeneous cosmology with compact spatial section is one, the spatial metrics of the Einstein- 
Rosen waves S'^ x R^, the Kantowski-Sachs universe with spatial metric di'^ = dz^ + a^dril, 
black hole spacetimes and the whole class of stationary axisymmetric metrics in relativistic 
astrophysics have EIRD=2. Similarly, the (imaginary time) finite temperature theory E? x S^ 
has EIRD three and the Kaluza-Klein theory on M^ x B^ with compact internal space B has 
EIRD four. 

It is important to recognize that this infrared dimensional reduction feature is not due directly 
to the topology of the underlying space, but rather the 'finite size' of the system. Spaces which 
are not direct product spaces but which can approach product spaces in the limit of extreme 
deformations or large perturbations have similar IR behavior. In the example of direct product 
spaces shown above, S'^ x 5^ can be shown to be the extreme deformation limit of S^\ The 
eigenfunctions of S^ which are the hyper spherical harmonics Ynem{x^G,(j)) can be expressed as 
a product of the spherical harmonics YnHmiJ^ ■, (t^) and the Gegenbauer polynomial Gni{x)- The 
eigenfunctions e*"-^ on S^ we saw above are obtained as limits of G'„(x) as the angular momentum 
of 5"^ is decoupled from the "radial" equations governing Gne- Another example given in [47J is 
the (static) Taub universe which is a deformed Einstein Universe. An eigenvalue analysis there 
provides a vivid illustration of this infrared dimensional reduction effect. To the extent that the 
large scale features of a quantum system shows up near criticality this reveals the surprising 
effect of its underlying spacetime structure. 

6. Summary 

In our first paper we explored macroscopic quantum phenomena (MQP) from a large N 
perspective, where N is the number of particles or components in a many-body quantum system. 
In this paper we explore MQP from the perspective of correlations and interactions which exist 
between components of the quantum system, how the coupling strength varies (runs) with 
energy or scale, especially how its infrared behavior at the critical point determines or affects 
the macroscopic features of a quantum system. We began with a review of the way how N, the 
number of particles in a quantum system, enters in the thermodynamic functions in quantum 
statistical mechanics. The (grand) partition function of a system at finite temperature (and 
chemical potential) in the (grand) canonical ensemble contains an explicit dependence of N, 
with proper factors accounting for whether the particles are distinguishable or identical. We 
note that the meaning of quantum in QSM refers specifically to spin-statistics (boson or fermions) 
and distinguishability, whereby classical and quantum are clearly divided. This restricted sense 
of quantumness in QSM is different from the other meanings or criteria of 'quantumness' which 
define MQP, i.e., quantum coherence, fiuctuations, correlation and entanglement associated with 
the full many-body wave function of a quantum system. 



We then explore the issue of how MQP can be measured or understood from the behavior 
of quantum correlations which exist in a quantum system and how the interaction strengths 
change with energy or scale, under ordinary situations and when the system is near its critical 
point. Noting that quantum correlation is further related to quantum fluctuations and quantum 
coherence, we used three examples from known results to illustrate these aspects: a) the existence 
of an H-theorem at the next-to-leading order large A^ expansion for an 0{N) quantum mechanical 
model, with entropy generation signifying the emergence of macroscopic (thermodynamic) 
properties, b) the nonequilibrium dynamics of A'^ atoms in an optical lattice under the large 
J\f (field components), 2PI and second order perturbative expansions, illustrating how N and 
J\f enter in these three aspects of quantum correlations, fluctuations and coupling strength, 
c) the infrared behavior of an interacting quantum system, we discuss the conditions where 
dimensional reduction shows up. The effective IR dimension is determined by the spectrum 
of the Laplacian of an interacting quantum system, in particular, whether a lowest eigenmode 
exists. On the same theme, we also discuss how the effective fleld theory concept bears on MQP: 
the running of the coupling parameters with energy or scale imparts a dynamical-dependent and 
an interaction-sensitive definition of 'macroscopia'. 

Quantum entanglement will be the theme of our third essay where we intend to focus on 
how this uniquely quantum feature shows up in quantum many-body systems and inquire to 
what extent it quantifies a quantum system as micro, meso versus macro. In this enquiry we 
will necessarily also touch on related subjects such as 'small systems' and emergent concepts in 
the new field known as quantum thermodynamics (ostensibly this is not the thermodynamics 
of quantum systems, treated in QSM textbooks, where quantum pertains to spin-statistics of 
particles). We can see why this is relevant to MQP in a simplified way: thermodynamics 
is a macroscopic theory valid when the volume V and number of constituents A^ both go to 
infinity while their ratio N/Y remains constant. Quantum dynamics is used for the description 
of micro-objects or small systems. As one increases the 'size' (number and volume, say) of a 
quantum system, in what parameter ranges will a thermodynamic description of such a many 
body quantum system become useful, at least begin to make sense? These inquires into the 
meaning of MQP will unavoidably lead us to confront some foundational issues of quantum and 
statistical mechanics. 

Acknowledgments 

We thank Esteban Calzetta for discussions at the initial stage of this work. BLH enjoyed the 
warm hospitality of the National Center for Theoretical Sciences (South) and the Department of 
Physics of National Cheng-Kung University, Tainan, Taiwan while part of this work was done. 
This work is supported by the National Science Council of Taiwan under grant NSC 97-2112- 
M-006-004-MY3, the US NSF grants PHY-0801368 to the University of Maryland. Two of the 
authors spent some time at the KITPC program on Quantum Open Systems supported by the 
Project of Knowledge Innovation Program (PKIP) of Chinese Academy of Sciences, Grant No. 
KJCX2.YW.W10. 



[1] 
[2] 
[3] 
[4] 

[5 

[6- 



Chou C H, Hu B L and Suba§iY 2011 Preprint quant-ph/1106.0556vl 

Huang K 1987 Statistical Mechanics 2nd ed. (Wiley) 

Zurek W H 1991 Phys. Today 44 36 

Guilini D et al 2003 Decoherence and the Appearance of a Classical World in Quantum Theory 2nd ed. 

(Berlin: Springer- Verlag) 
Pathria R K and Beale P D 1996 Statistical Mechanics 2nd ed. (Butterworth-Heinemann) 
Calzetta E, Ho Kwan-yuet and Hu B L 2010 Vortex Formation in Two-Dimensional Bose Gas, J. Phys. B: 

At. Mol. Opt. Phys. 43, 095004. 

Hu B L 2009 J. Phys. Conf. Ser. 174 012015 Preprint gr-qc/0903087.^ 

Calzetta E, Hu B L and Mazzitelli F D 2001 Physics Report 352 459-520 Preprtnt jhep-th/01 0219"9| 

Hu B L 1991 Coarse- Graining and Backreaction in Inflationary and Minisuperspace Cosmology Invited 



Lectures at the Seventh International Latin- American Symposium on General Relativity (SILARG VII) 

Proceeding appeared as Relativity and Gravitation: Classical and Quantum J. D'Olivio et al (eds) 

(Singapore: World Scientific) 
[10] Hu B L 2005 The Universe as an Ultimate Macroscopic Quantum Phenomenon, Invited talk at the Quantum 

Physics of Nature (QUPON) Conference Vienna., Austria 
[11] Jacobson T A 1995 Introduction to Black Hole Microscopy Utrecht THU-95/23, Preprint |hep-t h/951002"6] 
[12] Calzetta E and Hu B L 1993 Decoherence of correlation histories Directions in General Relativity vol. 2 Hu 

B L and Jacobson T A (eds) (Cambridge: Cambridge Univ. Press) Preprmf |hep-th/9501040| 
[13] Verstraete F, Popp M and Cirac J I 2004 Phys. Rev. Lett. 92 027901 
[14] Calzetta E and Hu B L 1999 Phys. Rev. D 61 025012 
[15] Calzetta E and Hu B L 1995 Correlations, decoherence, dissipation and noise in quantum field theory Heat 

Kernel Techniques and Quantum Gravity vol. 4 of Discourses in Mathematics and Its Applications ed 

FulUng S A (A&M Univ. Press, College Station, TX) Preprint |hep-th/9501040| 
[16] Calzetta E and Hu B L 2000 Phys. Rev. D 61 025012 
[17] Calzetta E A and Hu B L 2003 Phys. Rev. D 68 065027 

[18] Mihaila B, Athan T, Cooper F, Dawson J and Habib S 2000 Phys. Rev. D 62 125015 
[19] Mihaila B, Dawson J and Cooper F 2001 Phys. Rev. D 63 096003 
[20] Balescu R 1975 Equilibrium and Nonequilibrium Statistical Mechanics (New York: John Wiley & Sons) 

Balescu R 1997 Statistical Dynamics, Matter out of Equilibrium (London: Imperial College Press) 
[21] Baym G and Kadanoff L 1962 Quantum Statistical Mechanics (New York: Benjamin) 
[22] Calzetta E and Hu B L 1988 Phys. Rev. D 37 2878 
[23] Schwinger J 1960 Field Theory Methods in non-field-theory contexts Brandeis Summer Institute in Theoretical 

Physics. Keldysh L V 1964 Zh. Eksp. Teor. Fiz. 47 1515 (1965 Sov. Phys. JETP 20 1018) 
[24] Cornwall J M, Jackiw R and Tomboulis E 1974 Phys. Rev. D 10 2428 
[25] Hu B L 1989 Physica A 158 399 
[26] Hu B L 1994 Quantum statistical fields in gravitation and cosmology Proc. Third International Workshop 

on Thermal Field Theory and Applications Kobes R and Kunstatter (eds) (Singapore: World Scientific) 

Preprint gr-qc/9403061 
[27] Hu B L and Sinha S 1995 Phys. Rev. D 51 1587 
[28] Kac M and Logan J 1979 Fluctuations Fluctuation Phenomena eds MontroU E and Lebowitz J (New York: 

Elsevier) 
[29] Spohn H 1991 Large Scale Dynamics of Interacting Particles (Berlin: Springer- Verlag) 
[30] Habib S 1990 Phys. Rev. D 42 2566 

[31] Aarts G, Ahrensmeier D, Baier R, Berges J and Serreau J 2002 Phys. Rev. D 66 045008 
[32] Hu B L and Pavon D 1986 Phys. Lett. B 180 329 
[33] Habib S, Kluger Y, Mottola E and Paz J 1996 Phys. Rev. Lett. 76 4660 



[34] Habib S 2004 Gaussian dynamics is classical dynamics Preprint quant-ph/0406011 



[35] Calzetta E and Hu B L 2008 Nonequilibrium Quantum Field Theory (Cambridge: Cambridge University 

Press) 
[36] Fisher M P A, Weichman P B, Grinstein G and Fisher D S 1989 Phys. Rev. B 40 546 

Oosten D van, Straten P van der and Stoof H T C 2001 Phys. Rev. A 63 053601 
[37] Rey A M, Hu B L, Calzetta E, Roura A and Clark C W 2004 Phys. Rev. A 69 033610 
[38] Morgan S A 2000 J. Phys B 33 3847 
[39] Hutchinson D A W ef al. 2000 J. Phys B 57 3825 
[40] Hohenberg P and Martin P 1965 Ann. Phys. (N. Y.) 34 291 
[41] Shi H and Grifiin A 1998 Phys. Rep 304 1 

Griffin A 1996 Phys. Rev. B 53 9341 
[42] Giorgini S 1998 Phys. Rev. A 57 2949 

Giorgini S 2000 Phys. Rev. A 61 063615 
[43] Proukakis N P and Burnett K 1995 J. Res. Natl. Inst. Stand. Technol. 101 457 

Proukakis N P, Burnett K and Stoof H T C 1998 Phys. Rev. A 57 1230 

Kohler T and Burnett K 2002 Phys. Rev. A 65 033601 
[44] Calzetta E, Hu B L and Ramsey S A 2000 Phys. Rev. D 61 125013 
[45] Wright E M, Wong T, CoUett M J, Tan S M and Walls D F 1997 Phys. Rev. A 56 591 
[46] Smerzi A and Raghavan S 2000 Phys. Rev. A 61 063601 
[47] Hu B L and O'Conner D J 1987 Phys. Rev. D 36 1701 
[48] Hu B L 1987 Dynamical Finite Size Effect, Inflationary Cosmology and Thermal Particle Production, invited 

talk given at the CAP-NSERC Summer Institute in Theoretical Physics, Edmonton, Canada. Proceedings 

edited by K. Khanna, G. Kunstatter and H. Umezawa (Singapore: World Scientific Publishing Co. 1988). 



